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Abstract 



ON The effect of non-locality in the NN interaction with an off-energy shell character 

has been studied in the past in relation with the possibility that some models could 
be approximately phase-shifts equivalent. This work is extended to a non-locality 
implying terms that involve an anticommutator with the operator p 2 . It includes 
both scalar and tensor components. The most recent "high accuracy" models are 
considered in the analysis. After studying the deuteron wave functions, electro- 
magnetic properties of various models are compared with the idea that these ones 
differ by their non-locality but are equivalent up to a unitary transformation. It is 
found that the extra non-local tensor interaction considered in this work tends to 
re-enforce the role of the term considered in previous works, allowing one to explain 
almost completely the difference in the deuteron D-state probabilities evidenced by 
the comparison of the Bonn-QB and Paris models for instance. Conclusions for the 
effect of the non-local scalar interaction are not so clear. In many cases, it was found 
that these terms could explain part of the differences that the comparison of predic- 
tions for various models evidences but cases where they could not were also found. 
Some of these last ones have been analyzed in order to pointing out the origin of 
the failure. 
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1 Introduction 



A major aim of electron scattering experiment on few-body systems is to val- 
idate the ingredients entering the description of the process. These include 
the NN interaction models and the associated currents. With this respect, 
the recent measurements of the deuteron form factors at JLab [1,2] represent 
a valuable information. The interpretation of these data is not an easy one 
however. There are uncertainties on the nucleon form factors (see Ref. [3] for 
a recent review). The effect of relativistic corrections, which begin to show 
up in the explored Q 2 range, is approach dependent for a part which is not 
negligeable [4]. Some genuine effects, due to two-pion exchange, to an under- 
lying quark structure of the constituents, to the field-theory foundation of the 
interaction are expected to show up but they are likely to be mixed with 
other ones. Moreover, within a given scheme (corresponding to some equation 
for instance), NN interaction models, seemingly different, may represent the 
same physics. This involves different non-localities as expected from phase- 
shift equivalent potentials obtained by unitary transformations. 

As is well known, within a non-relativistic approach, it is always possible to 
describe the NN interaction below the meson-production threshold by a local 
potential for each partial wave separately. However, theoretically, the interac- 
tion is expected to be non-local. This observation motivated a recent study 
on the role of the non- locality, essentially on a phenomenological basis [5]. 
The non-locality in interaction models has many sources and often results 
from accounting for a theoretical ingredient, explaining that it is not always 
emphasized by the authors. Part of it originates from relativistic corrections 
(y / m/e factors, structure of Dirac spinors). It has been considered in Bonn 
models [6] and is sometimes approximated by linearly p 2 /m 2 dependent terms 
which makes them tractable in configuration space (Paris [7], Bonn-R [6] and 
Nijmegen [8] models). Another part comes from the coupling of the dominant 
NN channel to other ones, AA or nNN components. In these cases, the orig- 
inal time non-locality is transformed into a spatial one. This non-locality is 
not explicitly present in the models (except of course for the energy-dependent 
ones, full-Bonn for instance) but may be accounted for in the fitted param- 
eters. Finally, the compositeness of the nucleon or the exchanged meson can 
also produce some non-locality. It is explicit in quark based calculations (see 
Ref. [9] for a recent work together with earlier references) but, most often, 
as in the previous case, it is hidden in the parametrization of the potential 
model. 

While different non-localities are present in various models, it is not clear 
whether it helps in describing the NN interaction. These local and non-local 
models may be as many different models, in which case discriminating be- 
tween them is important. One cannot exclude however that they represent 
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the same physics for some part. In this case, there is not so much to study at 
the NN level and, most probably, one should look at other observables. In the 
example of form factors, a preferred model would require less corrections from 
two-body currents beside the genuine ones. The possibility to have different 
models representing the same physics was discussed in the past by Friar [10,11]. 
For a long time, it was just a curiosity till it was realized that more accurate 
models, namely Paris [7] and Bonn-QB [6], could realize to some extent this 
conjecture [12,13]. Terms under consideration have an off-shell character and 
mainly concern the tensor component of the interaction. When these models 
were used with the corresponding currents, it turned out that the discrepan- 
cies between the predictions for electromagnetic observables tend to decrease. 
Some was left however at the same time that some difference in the deuteron 
D-state probability was unexplained. It was mentioned in Refs. [12,13] that 
some role could be played by terms which have the form of an anticommutator 
of the kinetic energy with a tensor operator, p 2 W T + W T p 2 - Interestingly, Von 
Geramb et al. [14] , using inverse scattering methods, were able to derive local 
equivalent potentials for the two Paris and Bonn-QB models, evidencing very 
close deuteron D-state probabilities, suggesting that the above term could also 
be important in making the predictions of the two models closer to each other. 
These results largely motivate the present paper. 

Besides non-local contributions of the type, i(p 2 U — U p 2 ), which have been 
studied in a previous paper and later on by de Forest [15], we therefore consider 
non-local contributions of the type, p 2 W + W p 2 . As in the other case, W can 
have a scalar as well as a tensor component. The first one may originate 
from overall relativistic \jvaje normalization factors [16], mostly ignored in 
non-relativistic approaches. Due to the cancellation between the attractive 
and the repulsive components of the force due to a (2n) and uj exchanges, 
its effect is not expected to be very large. A larger effect may be due to 
the different Lorentz structures of these couplings, implying that the p 2 /m 2 
corrections to each of them add together instead of cancelling. The effect, 
which in relativistic mean-field approaches tends to decrease the a exchange 
contribution has a net repulsion. It is part of many models [7,6,8] and we 
believe that it is theoretically relevant. Its size is not well known however 
and, as far as we can see, there was no systematic attempt until now to look 
for an optimized strength. Perhaps, a reason for this is that such a term can 
always be transformed away by a unitary transformation as we recalled. If it 
has some relevance, one should expect that a smaller number of parameters be 
necessary to fit the same body of data. An indication in this sense is given by 
the comparison of the two Nijmegen models, Nijl and Nij2 [8]. Approximately 
based on the same philosophy (vertex form factors, ...), but differing by a 
linearly p 2 dependent term, they reproduce equally well NN scattering data 
(x 2 per datum = 1.03). The first one, which is non-local, requires however 
less parameters (41 instead of 47). The non-local component with a tensor 
character is primarily due to the pion exchange. This is a well established one 
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as far as the long-range force is concerned. 

When studying the effect of non-local terms of the form, p 2 W + Wp 2 , we 
will consider a transformation up to first order in the interaction. Contrary to 
off-shell effects generated by the terms, i(p 2 U — Up 2 ), the determination of 
the unitary transformation at the lowest order is not straightforward. It re- 
quires to solve an integral equation, which can only be performed numerically. 
When this is done, it is possible to calculate the effect on observables and 
see whether such terms can explain discrepancies between models evidencing 
different non-localities. This concerns the comparison between models. For 
our purpose, it is not necessary to look at the corrections to the interaction 
model which result from eliminating some non-locality. We assume that the 
fit of the model to scattering data has accounted for them, at least to a good 
approximation. We nevertheless notice that such corrections were explicitly 
considered in the past (see Ref. [17] for a recent work). Ultimately, one is in- 
terested in a comparison with experiment. With this last respect, the previous 
comparison may be useful to discriminate between effects that are simply due 
to the mathematical representation of the model and those due to genuine as- 
pects such as a two-pion-exchange contribution [6,7], the choice of the strong 
nNN coupling [18,19], ... 

The plan of the paper will be as follows. In the second section, we present the 
non-local terms considered in this work and show how they can be transformed 
away by a unitary transformation at the lowest order. This involves two terms 
with a scalar and a tensor character. The effect on wave functions is discussed 
in the third section for the deuteron case, including both its S and D state 
components. The fourth section is devoted to the effect on the deuteron form 
factors and the structure functions, A{Q 2 ) and B{Q 2 ) as well as the tensor 
polarization, T 2 o(Q 2 ). This part incorporates effects with an off-shell character 
described in an earlier work. Comparison of predictions for different interaction 
models is made when correcting them for a difference in their non-locality. A 
few cases where the discrepancy between predictions is not explained, and 
sometimes enhanced, are examined. A discussion and a conclusion are given 
in the fifth section. 



2 Formalism 

Non-locality can be produced by various contributions to the interaction. That 
one we are considering here originates from an expansion of the spinors en- 
tering the meson-exchange interaction. The part, which is by far the best 
identified, is due to the pion exchange. Its center of mass expression, up to 
first order terms in p 2 /M 2 , is given by: 
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The above equation may be written as: 



(P" - P 2 ) ■ (1) 



V = Vg + V? + {£ W£} + W?} + , i IT], (2) 



It can be extended to the exchange of other mesons: 

2 2 2 

V = V S + V T + Ws} + {fj, *M + < £/]. (3) 



In the above expression, Ws and WV can be expanded in such a way to exhibit 
their spin-isospin structure. In configuration space, they would read: 

s,t 

W T = J2 W T(r)S 12 (r)P lti , (4) 
t 

where the operators P s ^ represent the projectors on different spin and isospin 
channels. The spin tensor operator, shortly written S12, is defined as Si 2 (r) = 
a[ ■ f &2 • r — (f! • (T2/3. Together with this definition, we write the deuteron 
wave function as \ip(r) > oc [u(r) + 3w(r)/y / 8 S 12 (f)] \ 3 Si >. Working with 
states of given spin and isospin, these projectors can be omitted as well as the 
corresponding indices. The spin-tensor operator cannot be however. Without 
loss of generality, we can therefore make the replacements: 

W s <-> W s ( = W s {r) ), W T <-> W T S 12 ( = W T (r) S 12 (r) ). (5) 



The two first terms in Eq. (3), V5 and Vr, represent the standard scalar and 
tensor parts of the non-relativistic interaction. These terms are local and are 
present in all interaction models. They differ by the precise value of the cou- 
pling constant, g n NN, and the expression of the hadronic form factor. This one 
is taken as 1 like in the Reid models or assumes a more realistic form like in 
most other models. The potentials that only retain the two first terms at the 
r.h.s. of Eqs. (2, 3) are considered as local ones. They in particular include the 
so-called RSC[20], Reid93, Nij2 [8] and Argonne V18 [21] models that will be 
evoked in this work. They are denoted here as type-I models. 
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The so-called Paris [7], Nijl and Nij93 [8] models contain a minimal non- 
locality given by a term like the third one in Eq. (3). Partly justified by 
its appearance at the first order in p 2 /M 2 , the limitation to this term is also 
motivated by the fact that it can be dealt with relatively easily in configuration 
space. These models, like the type-I ones, are generally built in configuration 
space. They are denoted as type-II models. Interestingly, these models tend 
to show a suppression of the deuteron S-wave at short distances smaller than 
local models corresponding to a similar description of the short range part of 
the interaction. An open question is to know whether this can be related to 
the non-local term proportional to W s in Eq. (3). 

A last class of models (denoted type-Ill) is represented by the Bonn ones, QA, 
QB, QC, CD [6]. Contrary to the models discussed above, they consider the 
full structure of the Dirac spinors describing on mass-shell nucleons. Due to 
this, these models contain all kinds of non-locality present in Eq. (2) for the 
simplest pion-exchange case and in Eq. (3) for the most general case. They in 
particular involve terms that have an off-shell character and were studied in 
earlier works [12,13,15]. The dominant contribution was provided by the tensor 
part of U n which could explain 2/3 of the difference in the deuteron D-state 
propabilities evidenced by the comparison of calculations performed with the 
Paris and Bonn-QB models. An open question here, indirectly answered by 
von Geramb [14], is to know whether the fourth term in Eq. (3), oc Wt, can 
explain the remaining part. 

In comparing models, it is important that they are physically equivalent with 
respect to NN scattering properties. The so-called "high accuracy" models, 
Nijl, Nij2, Argonne V18 and Bonn-CD should therefore retain our attention. 
Concerning effects in relation with the tensor part of the interaction, we will 
also include the Paris and Bonn-QB models which were used in a previous 
work. These two models provide values for the mixing angle, e±, which are 
relatively close to each other on the one hand, and close to the experimental 
ones on the other. 

While the good quality of the fit to NN scattering data is essential to avoid 
biases in the analysis we are going to perform, the failure to establish some 
relationship between the various in principle phase-shift equivalent models can 
have different origin. Beyond higher order corrections, it may reveal a genuine 
physical effect. Although it can be largely accounted for in the model fitting, 
the two-pion exchange contribution explicitly incorporated in the Paris model 
could be an example of such an effect. It may also reveal that some ingredi- 
ent in a model is, on the contrary, physically irrelevant. Moreover, but not 
independently, it may suggest that models are not completely determined by 
scattering data in the range where they are available. Looking for possible fail- 
ures can thus provide information on the underlying physics, which, after all, 
is the ultimate goal for studying the NN interaction properties. Anticipating 
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on the following, this aspect will mainly involve its scalar part. In the next 
subsections, we show how non-local terms present in Eq. (3) can be trans- 
formed away at the lowest order in p 2 /M 2 . This is done successively for the 
scalar and tensor components of the interaction. We will concentrate on the 
terms that have an anticommutator character, {p 2 /M, Ws} and {p 2 /M, Wt} 
in Eq. (3). The other terms of the commutator type, [p 2 /M, iU], which are 
relatively easy to deal with, were discussed elsewhere [12,13,15]. 

For definiteness, we denote by H and H' the two Hamiltonians that, in the 
ideal case, would be related to each other by a unitary transformation: 

H' = e s He s+ , (6) 



where S = —S + . The Hamiltonian, H, will be that one that contains a non- 
local term, Vnl which we want to transform away by making a unitary trans- 
formation and H' will be the resulting Hamiltonian. The two Hamiltonians 
therefore assume the following expressions: 



2 

H= M +V + Vnl - (7) 

ff = £ + V, (8) 

When the non-local term, Vnl, is transformed away to make it local, the 
interaction V acquires an extra contribution A V, defined as: 

A V = V - V. (9) 



Starting from a two-body interaction, A V will generally involve two-body but 
also many-body contributions. As to the solutions of the Hamiltonians, H and 
H', they are denoted \P and 



2.1 Non-locality in the scalar part of the interaction models 



In this subsection, we consider two models that differ by a non-local term 
proportional to Ws- The non-local term, Vnl, therefore assumes the following 
expression: 

V N L = {^,W s (r)}. (10) 

The physical origin of the last term in Eq. (10) has been briefly mentioned in 
the introduction. The first candidate for such a contribution is provided by the 
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examination of the uj and a (2 7r in a S-state) exchanges. The cancellation that 
occurs for these contributions at the zeroth order in p/M is absent at the order 
(p/M) 2 , hence a possible sizeable term. In comparison, the contribution due to 
the pion exchange in Eq. (2) contains extra factors of the order (p/M) 2 . Due 



to the above cancellation, minimal relativity factors ym/e, whose effect was 
studied in a first step by Elster et al. [16], may only be a part of the dominant 
contribution. Non-locality could also arise from eliminating degrees of freedom 
involving mesons in flight or baryons excitations. It however appears that 
the corresponding spatial and time non-localities, individually large, tend to 
cancel each other. This probably explains the success of a local potential in 
reproducing the spectrum of the Bethe-Salpeter equation [22]. 

Equation (6) together with the expression of Vjvl given by Eq. (10) cannot 
be solved easily but taking into account that the term we want to remove 
is of order (p/M) 2 , we can solve it order by order by making an expansion 
of the unitary transformation that we denote e Ss . At the lowest order in the 
interaction, we have thus to fulfill the following equation: 

{£,W 8 } + [S S £]=AV§, (11) 



where A Vg must be local. This is achieved by choosing S s as: 

Ss = i f fr V (r) + V (r) r-p s ^ 

where the function V satisfies the following first-order differential equation: 
V (r) + f-VV (r) = W s (r) . (13) 

The solution of this equation is given by the following integral: 

1 °° 

V (r) = -- J W s (r')dr'. (14) 



The integral can be calculated numerically provided that the integrand is not 
too singular. The singular character of V due to the presence of the front 
factor 1/r is partly removed by the presence of the factor, r, in the expression 
of Ss, Eq. (12). It is likely that this singular character is the consequence of 
limiting ourselves to the lowest order terms in both expanding e Ss and in the 
interaction. It should be removed in a complete calculation but this supposes 
a fully non-local unitary transformation. 
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While the non-local term in Eq. (3) is removed at the lowest order, an extra 
local contribution to the interaction with same order in the coupling con- 
stant is generated. Denoted AVg in Eq. (11), its expression is given in the 
appendix. Let's precise here that its locality holds for each partial wave sep- 
arately, as expected from phase-shift equivalent models. On the other hand, 
it is somewhat more singular than the interaction models we started from or 
contains terms with extra powers of 1/r. By fitting the model parameters to 
the scattering NN data, such terms may be partly accounted for but we do 
not expect to completely explain the discrepancy between different models, 
the new terms having a functional form that is generally discarded in building 
models. Finally, an important difference is to be noticed with the treatment 
of the non-local term of the commutator type, [p 2 /M, iU], made elsewhere 
[12,13]. The correction, AV, contains in the present case a contribution of the 
first order in the coupling g 2 , AV°, while it was of the second order for this 
other term. 

For practical applications, we will consider the two potential models, Nij 1 and 
Nij2, which have been built with approximately the same ingredients, except 
evidently for the presence of a non-local term in the first one, which has the 
form of the term considered in this work. We also considered other couples of 
models, Paris and RSC for instance. Their probing strength as to their pos- 
sible unitary equivalence was not found very strong. Instead, the comparison 
pointed out to a discrepancy having its source in the different values for the 
quantities, A s and a t (and possibly the ratio r\ = A D /A S ) obtained with these 
models. Some correlation between A s and a t is well known but a correlation to 
Ad/As, which appears to be supported by the examination of various models, 
is not so. Characterizing the asymptotic behavior, the difference evidenced by 
the comparison of these values that are observable ones, cannot be removed 
by the unitary transformation considered here. 



2.2 Non-locality in the tensor part of the interaction models. 



In this subsection, we consider two models that differ by a non-local tensor 
term of the anticommutator type: 

V NL = {^,W T (r)S l2 (r)}- (15) 



The non-local tensor term of the commutator type has been discussed at length 
in earlier works [12,13]. 

To eliminate the non-local term under consideration, we proceed as in the 
previous subsection, with the same approximations (expansion of the uni- 



9 



tary transformation, now denoted e T , up to first order in St and first-order 
terms in the interaction retained). The quantity St has to be chosen in such 
a way that its commutator with the kinetic energy term p 2 /M cancels the 
term {p 2 /M, W T S 12 } without generating the same terms as those we want 
to remove. This implies to solve an equation similar to Eq. (11) for the scalar 
case: 

2 2 

{j^,W T S 12 } + [S T ,^]=AV^, (16) 

where A has to be local. This supposes to take two terms in S T and is 
achieved by: 

S T = ^(p-r S l2 (f) Vi(r) 

+ (a[-p a 2 -f + a 2 -p ux-r - | a r a 2 p-rj V 2 {r) + h.c^j , (17) 
The functions, Vx(r) and V 2 (r), have to satisfy the following relations: 



r(Vl(r) + V 2 \r)) = W T (r), (18) 

Vx(r) + V 2 (r) + (rV 2 (r))' = 0. (19) 

Despite it is more complicated than Eq. (13), the above system of equations 
can be formally solved with the result: 



oo oo . 



r r i 

1 7w T (r') , 

Vx{r) = -V 2 {r)--j-?ydr>. (21) 

r 

Numerical values for V\{r) and V 2 (r) are then obtained by performing suc- 
cessively the integrations in Eqs. (20) and (21). Concerning the singularities 
entering these expressions, the same observations as for Vb(r) can be made 
here. 

Similarly to the scalar case, the removal of the non-local tensor term of the 
anticommutator type generates an extra contribution to the interaction, AV^ , 
with a local character. Its expression is given in the appendix. Again, it evi- 
dences an analytic structure different from what is usually included in poten- 
tial models. A large part of it is certainly accounted for by fitting its parame- 
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ters to scattering data. Whether the remaining part can show up is probably 
difficult to answer in view of other uncertainties. 



Applications will be made for the couple of the Paris and Bonn-QB models that 
we considered in a previous work for part of the non-local tensor component 
with a commutator type, assuming implicitly that non-local scalar terms of 
the anticommutator type are more or less the same. We also considered the 
couple of models Argonne V18 and Bonn-CD. In this case, a non-local scalar 
term is also included, slightly obscuring the analysis. 



3 Comparison of wave functions 



In this section, we compare wave functions obtained from various models differ- 
ing by their non-locality. In principle, and as far as the models are phase-shift 
equivalent, it should be possible to relate the corresponding wave functions 
according to the equality: 

= e s (22) 



where \& and are solutions of Schrodinger equations with Hamiltonians 
respectively given by Eqs. (7) and (8). 

Consistently with our approximation retaining the first term in the expansion 
of the unitary transformation e s made in its determination, we assume that 
the correction to \1/ produced by the application of this transformation is given 
by 

= S *. (23) 



In practice, we will calculate this quantity, add it to \I/ and compare the result 
to the wave function obtained independently in a different model, W, which 
it could be compared with. Expressions of A\I> in the case of a non-local term 
with an off-shell character, oc iU, have been given in Refs. [12,13]. We here 
give these expressions for the other non-local terms, oc Ws and oc Wt- 

Correction to the wave function from transforming away the non-local term 
{p 2 /M,W s } 

For the first case, one gets separately for the S and D wave components of 
the deuteron wave function: 



Au = ^W s u + V ru', (24) 
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Aw = ^W s w + V rw'. (25) 

Using the expression of Vq in term of Ws, Eq. (13), one can easily check that 
the correction to the norm, always at the first order in S, is zero separately 
for the S and D waves: 

oo oo 

J dr u Au — J dr w Aw = 0. (26) 



This result is consistent with the unitary character of the transformation and 
the fact it does not mix the two waves. 

Correction to the wave function from transforming away the non-local term 
{p 2 /M, W T S 12 } 

Expressions for the case of the term oc Wt are more complicated: 



Am = ^ ((V 1 + 5V 2 + ^W T )w + (V 1 + 2V 2 )(rw' -w)^j, (27) 

A«; = y (Vx - V 2 + X - W T ) u + (V x + 2 V 2 ){ru' - «)) 

~ ({V 1 + 2V 2 + ^W T )w + (V 1 + 2V 2 )(rw' -w)j . (28) 



As previously, using the above expressions, one can check that the correction 
to the norm is zero but the result now holds for the total sum of the S and D 
waves contributions: 

oo 

J dr (uAu + w Aw) = 0. (29) 



Comparison of wave functions for Nijl and Nij2 models 
In Figs. 1 and 2, we show the effect of the transformation discussed above in 
the case of a scalar non-local term, {p 2 /M, Ws}. This is done for the couple 
of models Nijl and Nij2 that precisely differ by such a term. As far as the 
unitary transformation is truncated to the first order term in S, applying 
the transformation e s on \1/ or the inverse transformation e~ s on is not 
equivalent. We therefore present results where the two cases are considered 
in the same figure. Beyond 0.8 fm, where corrections are small enough to be 
treated perturbatively, it is found that both the S and D waves become closer 
to the wave functions they should be compared with (see the insert). Some 
discrepancy remains at the largest distances. In the D-wave case, it partly 
reflects a difference in the asymptotic normalizations, Ad, that are not affected 
by the unitary transformation. In the S-wave case, it is not clear whether the 
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Fig. 1. Effect of a non-local scalar term, {p 2 /M, Ws}, on the deuteron S-waves: 
The effect is shown for the Nijl and Nij2 models. The corresponding deuteron 
S-waves, u(r), are represented by the continuous and dashed lines respectively. They 
should be compared with the transformed waves obtained from the other model, 
represented in the figure by dashed-dotted and dotted lines. The insert shows an 
enlarged part of the figure in the region where the non-locality effect shows up at 
large distances. For this calculation, Ws is taken from the Paris model. 



remaining difference around 1.4 fm is an indirect consequence of the previous 
one, of the choice of Ws (actually taken from the Paris model in a first step 
analysis) or simply of the model parametrization. It is noticed that the wave 
functions at these distances are likely to be sensitive to the medium range of 
the interaction, which itself involves scattering data at intermediate energies. 
The effect of a difference in reproducing these data, such as the 8% one for the 
ei mixing angle around 100 MeV for the two models considered here, cannot 
be removed, evidently, by the unitary transformation. 

At lower distances, examination of the figures, especially for the S-wave, evi- 
dences structures. They are due for a large part to the perturbative character 
of our approach. At the lowest distances, the effect of the transformation goes 
in the right direction. It therefore seems reasonable to attribute the smoother 
behavior of the S-wave for some models to their p 2 dependence. The effect is 
however too large. It is likely that a more accurate treatment would transfer 
part of the effect at larger distance, around 0.6 fm, where the present effect is 
too small or even has a sign opposite to what was expected. This transfer is 
not arbitrary however as the norm should remain unchanged. 
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Fig. 2. Effect of a non-local scalar term, {p 2 /M, Ws}, on the deuteron D-waves: 
The effect is shown for the Nijl and Nij2 models. Same as in Fig. 1, but for the 
deuteron D-waves, w(r). 




Nij1 

Nij1 transformed 
Nij2 

Nij2 transformed 




1.0 1.5 2.0 

r(fm) 

Fig. 3. Effect of a non-local scalar term, {p 2 /M, Ws}, on the deuteron S- and 
D-waves: The effect is shown for a different Ws in calculating the transformed 
S- and D- waves from the Nijl and Nij2 models. The corresponding curves are 
respectively given in the right and left parts; same as in Fig. 1 otherwise. For this 
calculation, Ws is taken from a meson-exchange model. 

The existence of the above structures may introduce a bias in the comparison 
of form factors and other observables that will be made in the next section. 
Actually, the effect is limited by the fact that in calculating the form factors, 
we only retain the linear corrections in order to preserve the norm, according 
to Eqs. (26) and (29). Moreover, in the integrand, this correction is multiplied 
by a wave function which tends to zero in the relevant domain. The actual 
effect on form factors would thus involve a factor g 2 r 3 /24 in the integrand 
and is consequently suppressed. It is roughly of the order q 2 rQ/9Q where r 
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r(fm) 

Fig. 4. Effect of a non-local tensor term together with off-shell contributions: The 
effect is considered for the S- waves of the Paris and Bonn-QB models. Curves show 
the quantity u{r)/r to better emphasize the short range behavior. The thin and thick 
continuous lines represent the corresponding quantity for these models. The trans- 
formed wave functions accounting for the non-local tensor term, {p 2 /M, Wt S12}, 
are given respectively by dotted and long-dashed lines (labelled NL). The wave func- 
tions accounting for the whole effect are given by the short-dashed and dash-dotted 
lines. These last ones should be compared with the wave function of the other model. 

is of the order of 0.4 fm. Furthermore, one cannot exclude cancellations with 
corrections at higher r. As will be seen in the next section, it is what occurs 
and the different structures should be considered consistently. 

Due to the theoretical uncertainty on Ws, we estimated the transformed wave 
functions with assuming that this quantity was produced by single meson 
exchanges (ir,p,u,a) with standard values for the parameters. Results are 
presented in Fig. 3. They evidence smaller structures than in Figs. 1 and 2, 
probably in relation with a smaller strength of the Ws term. By using in the 
following a term taken from the Paris model, it is likely that its strength is 
not underestimated. 

Comparison of wave functions for Bonn-QB and Paris models 
In Figs. 4 and 5, we show the effect of the transformation discussed above 
in the case of a tensor non-local term, {p 2 /M, W T S 12 }- For this calculation, 
Wt S12 is taken as the sum of a 7r-exchange contribution appearing in Eq. (1), 
which provides the dominant contribution, and a p-exchange contribution, 
which has a cut-off role. Together with this non-local term, we consider the 
effect of the off-shell term studied in a previous work. A monopole cut-off 
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Fig. 5. Effect of a non-local tensor term together with off-shell contributions: The 
effect is considered for the D-waves of the Paris and Bonn-QB models. Curves 
represent the usual D-waves, w(r); same as in Fig. 4 otherwise. 



form factor, (A 2 - /i 2 )/(A 2 + k 2 ), with A = 1.05 GeV, is introduced in all 
calculations [13] so that to kill the short-range part of the interaction, for which 
the analysis developped in this work is more uncertain, both theoretically and 
mathematically. For the S-wave (Fig. 4), the transformed wave functions are 
relatively close to the one they should be compared with. The effect of the first 
term, {p 2 /M, W T S 12 }, which arises from the action on the D-wave (therefore 
different of the scalar case) goes in the right direction. Most of the effect 
comes from the term with an off-shell character. Notice that for this curve, 
we plotted u(r)/r, in order to better see the effects at short distances. Quite 
similar results hold for the D-wave shown in Fig. 5. It is immediate from these 
results that when considering the electromagnetic properties of the deuteron, 
the predictions of the different models will tend to be close to each other when 
the currents that are associated to each model are accounted for. Concerning 
the little structures appearing in the transformed wave functions, the remarks 
made previously for the comparison of Nijl and Nij2 wave functions could be 
repeated here. The amplitude of these structures is much less important here 
however. 
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4 Effect on deuteron form factors and electromagnetic observables 



Given a physical system and an underlying theory, it can be described by an 
infinity of interaction models, depending on the choice of degrees of freedom. 
This is quite similar to the choice of a frame or a basis in making an experiment 
or a calculation. This freedom is at the origin of the non-observable character 
of the wave function. Concerning the system itself, this is realized by using 
phase-shift equivalent potentials. As for the interaction with an external probe, 
described by an operator denoted O, the equivalence of the models implies a 
relation such as: 



In this expression, the operator O', which can contain in the minimal case one- 
body terms but also many-body ones, will here refer to the electromagnetic 
interaction. Using the expression relating ty' and Eq. (22), it is found that 
the above equality supposes that: 



In the case where O is a one-body operator, it is immediately found that the 
correction at the r.h.s. has in the minimal case a two-body character. In order 
to make a relevant comparison of predictions of observables from different 
representations, one has therefore to account for the contribution due to AO 
when calculations are performed. This is what we are doing in this section. The 
question of making a comparison to measurements, which in any case supposes 
to consider the contribution of other currents, is evoked in the conclusion. 

For the off-shell correction considered in Refs. [12,13], the commutator AO 
could be calculated explicitly, leading to two-body currents that, except for 
a factor, had for some part a structure identical to the pair-term contribu- 
tion to the charge operator. The expression for the magnetic operator could 
not be related to any known two-body operator. A similar situation holds 
for the transformation considered here for terms of the form {p 2 /M, Ws,t}- 
One thus gets the following contribution for the isoscalar charge operator in 
configuration space at the lowest order : 




(30) 



0' = e s Oe- s = + AO, 
AO = [S,0] + .... 



(31) 



[S,0] 




+ (<f 1 -q(J2-r + (j2-qa 1 -r--a 1 -a2q-r) V 2 (r)\. (32) 
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In practice, calculations can be performed indistinctively in two ways. The 
first one, in the spirit of equivalent models, consists in calculating the matrix 
element of O + A O with the wave function, and compare it to the matrix 
element calculated with O together with the wave function, \1>. In the second 
one, the first matrix element is replaced by the matrix element of O with the 
transformed wave function, e~ s While in this later approach conclusions 
immediately stem from the examination of transformed w.f. discussed in the 
previous section, the former approach is not without interest. Indeed, some 
of the corrections to charge, quadrupole and magnetic form factors (for the 
scalar part Wg) have a relatively simple expression: 



AF C (Q 2 ) = G S E (Q 2 ) J dr XJl {x) Uu 2 (r) + w 2 (r)) V (r) 
o ^ 

+^(«(r) w(r) - -L w 2 (r)) ( Vl {r) + 2 V 2 (r)) 



AF Q (Q 2 ) = ^ GUQ 2 ) J dr 



Q 2 

u 

x ^(Sj 2 (x) - atfi(aO) (u(r) w(r) - w 2 (r)) V (r) 
+V8j 2 (x) (-w 2 (r) +w 2 (r) + -^=u(r) w(r)) V 2 (r) 

(3j 2 (a;) -xj^x)) 

x (u 2 (r) + - w 2 (r) - V2 u(r) w(rfj (Vi(r) + 2 V 2 (r 

00 

AF M (Q 2 ) = 1 G S E (Q 2 ) J dr 3j 2 (x) w 2 (r) V (r) 



+G S M (Q 2 ) Jdr (xj.ix) (u 2 (r) - l -w 2 {r)) 
^ 

+ \ (3j 2 (z) - x 3l (x)) (w 2 (r) + V2 u(r) w(r)fj V (r) + ... (33) 

where x = Q r/2. The above expressions can be also obtained from those using 
the transformed wave functions by performing an integral by parts and taking 
into account that the term ....|o° is zero. Contrary to these ones however, the 
integrand is quite smooth and does not evidence structures at short distances. 
Far to suggest numerical inaccuracies or errors in the formalism, the structures 
that we discussed in the previous section have to be taken as they are and, most 
important, dealt with consistently. The dots at the last line of the expression 
for the magnetic form factor contain terms with the derivatives of the wave 
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functions. Unlike the charge and quadrupole form factors, these ones cannot 
be removed using an integration by parts. 

When comparing observables, we will proceed in two ways, depending on which 
model predictions are compared with. This is dictated by the first order ap- 
proximation made in expanding the unitary operator e s , which is required for 
consistency, as explained above, but prevents to keep the unitary character 
of the transformation. In assessing the two sets of results, it is likely that the 
truth should be in between. We successively consider the static observables, 
the form factors and finish the section by mentioning a few peculiar results. 



4-1 Static observables 
Effect of {p 2 /M,W s } 

When comparing predictions for models differing by a term {p 2 /M, W$}, we 
found very small effects for the static observables, r 2 D and Qo, and obviously 
none for \xd (in relation in this case with the absence of change of Pd)- The 
non-local term being short-range-behaved, we rather expect effects in the dy- 
namical observables (form factors) as it can be inferred from the discussion 
of the wave functions in Sect. 3. Concerning the models, Nijl and Nij2, we 
looked at more specifically, static observables are quite close to each other, 
leaving little room for possible effects. The difference in the values of r 2 Dl 
0.004 fm 2 , can be traced back to a difference in the asymptotic normaliza- 
tions A s for 75% (models roughly fulfill the relation Ar m ~ 1.9 AAg). What 
we calculated, ~ 0.001 fm 2 , is consistent with the unexplained part. For the 
quadrupole moment, we found a correction that has a sign opposite to what 
would be required by the comparison of the two models. Again, the major part 
of the difference in this case is likely to have another origin (see our discus- 
sion of the D wave around 1.4 fm in Fig. 2, Sect. 3). As for the the magnetic 
moment, it is unchanged since the normalizations of the S- and D-waves are 
separately unchanged. 

Effect o/{p 2 /M, W T S 12 } and other terms 

Effects of the tensor term, {p 2 /M, WtS 12}, are significantly more important 
than for the scalar one. In a previous work [13], we could explain 2/3 of 
the difference in the D-state probabilities for the Bonn-QB and Paris models 
(0.78%) by considering a non-local tensor term with an off-shell character, 
[p 2 /M, iU\. By adding the above term, we can now explain a difference of 
0.74% (0.14% for the new term alone). The deuteron D-wave probability enters 
in the magnetic moment according to the relation Afi£> = 3/2 (^5 — 1/2) APp. 
The difference for the two models is Afip is 0.0045 nm while we calculate 
0.0043 nm. Another observable sensitive to the tensor force is the quadrupole 
moment. Considering the quantity Qd/(A s A d ), that may be less sensitive 
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to differences in As or left unchanged by the unitary transformation we 
are dealing with, we found a difference of 0.19 fm 3 or 0.20 fm 3 depending on 
which model is corrected. The difference in the predictions is 0.23 fm 3 . Finally, 
we explain 60% of the difference in the squared charge radius evidenced by 
the comparison of the two models, Bonn-QB and Paris. This last observable 
is directly related to the slope of the deuteron structure function, A(Q 2 ). 

As far as static observables are concerned, it thus appears that they are very 
close to each other once the effects of their supposed unitary equivalence are 
taken into account. 



4-2 Comparison of form factors with models Nijl and Nij2, effect of the non- 
local scalar term, {p 2 /M, W$} 

The effect of the non-local scalar term, {p 2 /M, Ws}, on electromagnetic ob- 
servables is shown in three figures for the case of the Nijl and Nij2 potentials 
models, which differ by such a term. Here and below, we expect that if this 
term has some relevance in explaining the difference in the wave functions of 
these models, then predictions for the observables should be closer to each 
other when the appropriate currents are used. The ratio of these predictions, 
that we represent in the figures, should therefore approach 1. 




Fig. 6. Effect of a non-local term, {p 2 /M, Ws}, on charge form factors: Curves 
represent the ratio of the Nijl to the Nij2 charge form factors from Q 2 = to 
Q 2 = 15 fm -2 for the left part and from Q 2 = 25 to Q 2 = 50 fm -2 for the right 
part. The ratio of the uncorrected predictions is given by the continuous line. The 
ratios with the transformed Nijl and Nij2 w.f. (or including the associated two-body 
currents) are respectively given by the dashed and the dash-dotted lines. 

The effect on charge form factors is seen in Fig. 6 that has been split into 
two parts, due to the appearance of a zero in the range 15-25 fm -2 . While the 
effect was hardly significant for static observables, examination of the figure 
indicates that it is quite relevant at higher Q 2 and goes in the right direction. 
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Fig. 7. Effect of a non-local term, {p 2 /M, Ws}, on quadrupole and magnetic form 
factors: Curves represent the ratio of the Nijl to the Nij2 quadrupole and magnetic 
form factor from Q 2 = to Q 2 = 50 fm~ 2 for the first one (left part) and Q 2 = to 
Q 2 = 30 fm -2 for the second one (right part). They are defined as in Fig. 6. 
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Fig. 8. Effect of a non-local term, {p 2 /M, Ws}, on A(Q 2 ) and T 20 (Q 2 ) observables: 
Curves represent the ratio of the Nijl to the Nij2 A(Q 2 ) and T2o(Q 2 ) observables, 
respectively for the left and right parts. They are defined as in Fig. 6. 

Results do not depend much on whether two-body currents are added to one 
model or to the other. A similar conclusion holds for the quadrupole form 
factor as well as the magnetic one shown in Fig. 7. Notice that the last one is 
given up to 30 fmT 2 due to the presence of a zero at a slightly larger Q 2 . 



It is also instructive to look at quantities which are closer to what is measured, 
namely the deuteron structure function A{Q 2 ), and the tensor polarization at 
9 e = 70°, T 20 (Q 2 ). They are given in Fig. 8. Not surprisingly, they nicely 
confirm the above conclusions. 



Results quite similar to the above ones have been obtained when extended to 
other interaction models (Nijl replaced by Paris or Nij93 ; Nij2 replaced by 
Argonne VI 8 or Reid93). The agreement is especially good for those models 
that better fit the NN scattering data. Due to the various uncertainties al- 
ready mentioned, we will mainly retain that the currents associated with the 
removal of a term, {p 2 /M, Ws}, in the interaction can reasonably account for 
differences in model predictions, both in sign and in magnitude. Some differ- 
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ences remain however (and hopefully) unexplained. They will be discussed at 
the end of this section. 



4-3 Comparison of form factors with models Bonn-QB and Paris, effect of 
the non-local tensor term 
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Fig. 9. Effect of a non-local term, {p 2 /M, Wt S 12} together with off-shell effects 
for the charge form factor: Curves represent the ratio of the Bonn-QB to the Paris 
charge form factors from Q 2 = to Q 2 = 15 fm~ 2 for the left part and from Q 2 = 25 
to Q 2 = 50fm~ 2 for the right part. The ratio of the uncorrected predictions is given 
by the continuous line. The ratio when the Paris (Bonn-QB) w.f. is transformed 
for the first ( {p 2 /M, Wt S '12}) and second (off-shell) non-locality is given by the 
dotted and short-dashed (long-dashed and dash-dotted) lines. 
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Fig. 10. Effect of a non-local term, {p 2 /M, Wt Su} together with off-shell effects 
for the quadrupole and magnetic form factors: Curves represent the ratio of the 
Bonn-QB to the Paris quadrupole and magnetic form factors (left and right parts 
respectively). They are defined as in Fig. 9. 

The effect of a tensor term like {p 2 /M, W T S 12 } on electromagnetic observ- 
ables is shown in three figures for the case of the Bonn-QB and Paris models, 
which differ by such a term. These models also differ by terms that have an 
off-shell character and were studied in an earlier work. In presenting the var- 
ious effects, we first present the new one considered here and then add the 
older one. 
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Fig. 11. Effect of a non-local term, {p 2 /M, Wt S'12} together with off-shell effects 
for A{Q 2 ) and T2o(Q 2 ) observables: Curves represent the ratio of the Bonn-QB to 
the Paris A(Q 2 ) and T2q{Q 2 ) observables. They are defined as in Fig. 9. 
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Fig. 12. Effect of a non-local term, {p 2 /M, Wt S 12} together with off-shell effects 
for the deuteron electrodisintegration form factor near threshold: Curves represent 
the ratio of the Bonn-QB to the Paris form factor from Q 2 = to Q 2 = 11 fm -2 
for the left part and from Q 2 = 15 to Q 2 = 40 fm -2 for the right part. The ratio 
of the uncorrected predictions is given by the continuous line. The ratio when the 
Paris (Bonn-QB) w.f. is transformed for the first ( {p 2 /M, Wt S'12}) and second 
(off-shell) non-locality is given by the dotted and short-dashed (long-dashed and 
dash-dotted) lines. 

The effect of the non-local tensor term on the charge form factor is presented 
in Fig. 9. It explains a little part of the difference between the Bonn-QB and 
Paris model predictions. By far, the dominant effect is due to the off-shell 
effect associated to the n- and p-exchange tensor force. Quite similar results 
hold for the quadrupole and magnetic form factors shown in Fig. 10 and, as 
a consequence, for the structure function, A(Q 2 ), and the tensor polarization, 
T2o{Q 2 ), both given in Fig. 11. In these two cases, possible discrepancy is 
completely removed below 15 fm~ . There may be still some beyond in relation 
with that one observed for the quadrupole form factor in Fig. 9. We notice 
that in dealing with the non-local tensor force, we omitted in the Bonn-QB 
model a term of the form: 



cr r p' x p <7 2 -p' x p . 



(34) 
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This one appears at the fourth order in p/M and, moreover, arising from 
vector-meson exchange, has a short-range behavior. In configuration space, 
this term gives rise to a <j\ • L a 2 • L term, which is present in the Paris 
model considered in this work, non-local terms that are similar to those we 
discussed and another non-local one, proportional to a[-p d^-p [8]. This last 
one could be dealt with by using a transformation identical to that given in 
Eq. (17), but with different relations of V\ and V2 to the potential strength (see 
appendix). At first sight, this new term could contribute in the right direction. 
However, in discussing its contribution, we should also include higher order 
contributions due to eliminating non-local terms at the lowest order, [S, Vnl\- 
These induced terms have an extra factor 1/M 2 (for the tensor part) but a 
range quite similar to that of vector- meson exchange. 

Most results presented above concern an elastic transition. Some have been 
obtained for an inelastic transition, namely the deuteron electro-disintegration 
near threshold which involves a matrix element between the deuteron 3 S± (+ 3 Di) 
state and the 1 5'o scattering state. They are presented in Fig. 12. Examination 
of these results shows features quite similar to the ones observed in other cases. 
In principle, this process allows one to get insight on the effect of a non-local 
term in a channel, 1 S'o, different from the deuteron one. However, as most of 
the effect is related to the tensor force, the conclusion one can draw from this 
case is limited. 

As for the case of the non-local scalar term considered in the previous sub- 
section, we looked at other models (Bonn-QB replaced by Bonn-CD ; Paris 
replaced by Nijl or Nij93). Again, a better agreement is obtained with models 
that better fit NN scattering data. In the detail, a few points deserve to be 
discussed however. 



4-4 A few peculiarities 

The studies presented in this work were originally performed with the idea that 
various NN interaction models could be equivalent (up to a unitary transfor- 
mation). With this respect, it is interesting that we can explain a large part 
of the differences between models. The question now arises to determine the 
limitations, if any. By construction, the models are not identical. Besides a 
difference in their non-locality, discussed throughout this work, there may be 
differences in the choice of various ingredients: coupling constants, vertex form 
factors, exchanged mesons together with their mass, quality of the fit to NN 
scattering data, number of free parameters, etc.... After the non- locality effects 
have been removed, the comparison of model predictions can thus reveal more 
clearly features in relation with these ingredients. We here consider some of 
them. They involve a few percent effects that may not be significant at first 
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sight but become so when appearing repeatedly. They also mainly concern 
the low momentum transfers, the high momentum transfer range being too 
sensitive to a variety of ingredients to make firm statements. 
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Fig. 13. Comparison to results involving the Nijmegen 93 interaction model, for the 
charge form factor (left part) and the A(Q 2 ) structure function (right part). 

In Fig. 13 (left part), we show an example where an approximate agreement 
between two models turns into a discrepancy when the influence of the non- 
local interaction is accounted for. It concerns the charge form factor calculated 
with the Bonn-QB and Nij93 models in the range Q 2 = to Q 2 = 15 fm -2 . 
Being alone, this example, which is given to illustrate our purpose, does not 
allow one to draw conclusions on one of the two models under consideration. 
However, after incorporating the non-locality effect, it appears that the charge 
form factor for the Nij93 model evidences the same pattern as for the other 
models. This can be seen in the comparison of the A(Q 2 ) structure function 
calculated with the Nij93 model and with the set of models Nijl, Nij2, Reid93, 
Argonne VI 8 which, as already mentioned, provide very close results when 
non- locality effects are taken into account. The ratio of these quantities, shown 
in Fig. 13 (right part), evidences a difference in the slope or even a bump in 
the range Q 2 = to Q 2 = 15 fnT 2 . This tends to suggest that the Nij93 model 
involves specific features that are absent in the other models. It is tempting to 
make a correlation with the fact that this model predicts a smaller radius than 
the other ones, whose effect could not be seen clearly without removing effects 
of the non-locality. As to the origin of this smaller radius, it could be due to 
the poorer quality of the fit to scattering NN data (only 15 parameters). Not 
independently, it could also be due to the presence in this model of a pomeron 
exchange contribution which is ignored in all the other models. Its long range 
is likely to produce specific features that could not be easily accounted for 
with shorter range contributions to the potential. 

Another particular feature shows up when comparing the A{Q 2 ) structure 
function calculated with the Paris model to that one calculated with one of 
models Nijl, Nij2, Reid93, Argonne V18 and Bonn-CD and the appropriate 
corrections for non-locality effects. The ratio, which is shown in Fig. 14 (right 
part), evidences a slight bump in the range Q 2 = to Q 2 = 15 fm" 2 . A 



25 



1.04 



Nij1 




Fig. 14. Comparison to results involving the Paris interaction model, for the 
quadrupole form factor (left part) and the A(Q 2 ) structure function (right part). 

similar but less pronounced effect is seen with the Bonn-QB model. Again, 
the difference in the slope suggests to look at the deuteron radius predicted 
by the models. It is thus found that the radius for the Paris and Bonn-QB 
models, after accounting for non-locality effects, are slightly larger than for the 
other models. Before wondering about a possible origin for this difference in the 
radius, we emphasize another feature that may have some relation. It concerns 
the quadrupole form factor and more specifically the quadrupole moment. In 
many cases, we could explain differences in the quadrupole moment by the 
difference in the non-locality of the models. This is not the case for the Paris 
model (and to a lesser extent for the Bonn-QB model) as it can be seen in 
Fig. 14 (left part). The quadrupole moment for the Paris model is always 
larger than for the other models by an amount of 3-4%. There is no direct 
relation to the radius discussed above although a larger value will certainly 
contribute to produce a larger quadrupole moment. An explanation could be 
looked for in the value of the 7rNN coupling constant which is larger for the 
Paris (and Bonn-QB) model than for the other models which, for some of 
them, have been built on the prejudice of a value around g% NN /4ir ~ 13.5 
[18]. The value for the Paris model is closer to that one advocated by Ericson 
et al., g^ NN /4:7i ~ 14.11 [19]. A larger coupling implies a larger attraction 
at large distances and therefore a larger probability to find the neutron and 
the proton far apart in the deuteron. In the case of the Paris model, one can 
add that the explicit two-pion exchange contribution is likely to also produce 
a larger attraction at long distances in comparison with an approximation 
where this contribution is treated as an exchange of an effective a with mass 
around 500 MeV. 



5 Conclusion 



In this work, we have studied the effect of non-local terms on the description 
of the NN interaction. We especially looked at terms which involve the anti- 
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commutator with the operator p 2 , namely {p 2 /M, Ws,t}- The present work 
completes another one where terms of the commutator form [p 2 /M, iU], with 
an off-shell character, were considered. In dealing with these terms, we relied 
on an explicit form for the transformation allowing to get rid of them, which 
furthermore was treated at the first order in the interaction. In comparing 
models, this approach is probably not as good as an approach based on in- 
verse scattering methods. It nevertheless has the advantage to be based on the 
explicit identification of non-local terms in interaction models and to provide 
an insight on the currents that have to be associated with the models. As far 
as we can see, the method is accurate enough for our purpose, especially in 
view of its qualitative character for a part. On the other hand, the study has 
been extended to other interaction models, allowing one to draw conclusions 
that the comparison of two models only could not. 

For our study of the newly considered terms, we took the scalar part, Ws, 
from the Paris model while the tensor part, Wt, was taken from the sum of a 
contribution due to pion exchange, rather well determined, and a short range 
contribution. It is reminded that there is no a priori guarantee that these non- 
local terms can explain differences in interaction models which could have 
other sources and in some cases could not be explained (when they involve 
differences in observables such the asymptotic normalizations, As and Ad). 

We found that the scalar part, Ws, has quite small contributions to static 
quantities, such as the charge radius or the quadrupole moment. To evidence 
them, one has to get rid of contributions due to different asymptotic normal- 
izations for instance, which compete with the effects we are looking at or are 
even larger. Some sizeable effect however shows up for dynamical observables 
such as form factors. It is difficult to make a definitive statement on the size of 
the effect in view of various uncertainties, but its sign and its magnitude can 
reasonably explain differences between models differing by the above term. 

Concerning the tensor part, Wt, the contribution to static properties is more 
important. It, in particular, explains a non-negligeable part of the difference in 
the deuteron D-state probabilities between the Bonn models, which account 
for the full structure of Dirac spinors, and the other ones. When added to 
the effect of non-local terms with an off-shell character, which is larger, not 
much difference remains to be explained. For dynamical observables, the effect 
of the above term goes also in the right direction and accounts for a half of 
what could not be previously explained with the off-shell term. It would be 
interesting to see what would do an extra non-local tensor term omitted in 
our analysis. This one has a short range behavior but also a different structure 
which would require further theoretical elaboration. 

While a large part of differences evidenced by the comparison of models with 
various non-localities can be ascribed to this feature, we would not like to 



27 



give the impression that every difference can be explained that way. Inter- 
estingly, after removing the above effects, some differences between models, 
that were masked until now, more clearly appear. The correlation between 
the slope of the deuteron structure function, A(Q 2 ), and the deuteron radius, 
which is expected to hold in any case at very small Q 2 , is found to extend to 
higher Q 2 , up to 10 fm -2 . Knowing that there is for local models an empirical 
relation between the radius and the asymptotic normalization, As, the above 
observations raises the question of which model is the best with this respect, 
taking into account that this quantity is not changed in the transformations 
we performed. 

The examination of the quadrupole form factor, especially at Q 2 = 0, where it 
identifies to the quadrupole moment, Qr>, also shows interesting features. For 
this one, contributions involving the non- locality can amount to 1%, which 
is less than the difference between the upper value obtained with the Paris 
model (0.279 fm 2 ) and the lower ones (around 0.270-0.271 fm 2 ), obtained with 
the Nijl, Nij2, Reid93, Bonn-CD models. Curiously, there is no direct relation 
to the value of the ratio, i] = A^/As, which mainly determines this moment. 
After the effect of the non-locality is accounted for, which mainly affects values 
obtained with the Bonn models with respect to the other ones, we observe that 
the quadrupole moments are exactly in the same order as the 77 values and 
not far to be proportional to them. Accounting for the non-locality effects 
therefore restores some regularity that was lost when considering predictions 
of Bonn models. 

Throughout this paper, we ignored comparison with measurements, being es- 
sentially concerned by the comparison of theoretical predictions. We believe 
it is important in discriminating between models to separate effects that are 
due simply to a mathematical representation of the interaction, with more or 
less non-locality, and effects that have their origin in genuine aspects of the 
models. In making a comparison to measurements, one has first to determine 
among the various approximately phase-shift equivalent descriptions of the 
NN interaction which one is the more realistic, in the sense of relying on de- 
grees of freedom that are as little effective as possible. Then, one has to also 
account for contributions from two-body currents corresponding to this opti- 
mal model. Some work along these lines has recently been done with the aim 
to extract the neutron charge form factor from the knowledge of the deuteron 
quadrupole form factor [23]. 

We are very grateful to R. Machleidt, C. Elster and R. Schiavella for providing 
us with material that has been employed here or has been useful in analyzing 
our own results. 

APPENDIX 
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Induced contributions to the interaction, AV 

When the non-local terms of the form {p 2 /M, W} are removed from an in- 
teraction model to transform it into a local potential, local terms of the same 
order in the interaction are generated. Their expression is given here for the 
scalar and tensor parts. For the first one, there is a factor 2 change with what 
we give in Ref . [13] . The new convention now agrees with the definition of the 
non-local term in Ref. [7]. For the second one, we give an expression which 
evidences its tensor character. Due to relations between V and Ws, or V\, 
V 2 and W T for the scalar and tensor non-local terms respectively (see Eqs. 
13, and 18, 19), other expressions could be obtained. Those given below are 
among the simplest we could find: 

~S 12 (T)(W^+ 12^), (36) 

where 

S 12 (A) = i (a v A a 2 -A + a 2 -A a v A) - ^a v a 2 I 2 , (37) 
with A = L or A = f. 

As can be noticed, removing non-local terms introduces terms dependent on 
the total angular momentum, in agreement with the expectation that the 
phase-shift equivalent potential is local, wave by wave however. This is espe- 
cially transparent for the scalar part. It is also noticed that the global range 
is not affected, contrary to the removing of off-shell effects, but the potential 
is somewhat more singular due to the presence of derivatives of Vo, V±, V 2 , 
or Ws and Wt in the above expressions. Particular caution may be therefore 
required in making qualitative statements concerning the effect of these terms. 
For instance, in case where Ws would be given by a function, exp(— /i r), the 
term involving the second derivative in Eq. (35) would produce a contribution 
that is attractive at large distances. However, it also produces a repulsive S(r) 
term that will show up at very short distances. This one, often discarded, can 
give rise to the suppression of the S-state wave function that local models 
often evidence in this range. Although the singular character of the correc- 
tion to the interaction does not necessarily prevent one to perform meaningful 
calculations, we felt better to cut-off the short-range part of the non-local in- 
teraction. In this way, we perhaps underestimate some of the effects but those 
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we calculate are thus determined by the longer range part of the force, which 
is physically better known. 



Contributions to the static moments 

We here give the corrections to the static moments that can be obtained a 
direct calculation (Pd) or from Eqs. (33) for the form factors in the limit 
Q -> (A < r 2 >, AQ). 



r— OO / 

AP D = ^- J dr ((w T (r) -6V 2 (rj) u(r) w{r) 
o ^ 

+2r (y 1 (r)+2V 2 (rj) u\r) w(r))j, (38) 

i 7 

A < r 2 > s = 2j drr2 V °( r ) ( w2 ( r ) + w2 ( r )) > ( 39 ) 



A < r 2 > T = ^- J drr 2 (vi(r) + 2 V 2 (r) 
o 

x (•u(r) w(r) - -^=w 2 (r)), (40) 
AQ s D = y^J dr r 2 \/ (r) («(r) «;(r) - -j= w 2 (r)) , (41) 



2 7 ( 1 
AQ£ = -jdrr 2 i V 2 (r) (w 2 (r) - u 2 {r) + -= u{r) w(r)) 

o ^ ^ 

-l(y 1 (r) + 2V 2 (r)) (u 2 (r) + ^w 2 (r) - V2u(r) w(r))\. (42) 



About another non-local term 

We here consider a non-local term of the form: 



= i ar I ^ JL ^ w] (43) 



This one may arise from approximating the term given by Eq. (34) when 
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treated in configuration space (order 1/M 4 ). It could also come from the ex- 
change of a ai meson at the order 1/M 2 . It contains both a scalar part, which 
has a form identical to the one discussed in the main text of the paper and a 
tensor part. This last term can be transformed into a local one at the lowest 
order, using the same transformation as for the other non-local tensor term 
(17). The functions V\ and V 2 now fulfill the relation: 



Vl(r) + Vi(r) = 0, (44) 
V 1 (r) + V 2 {r) + (rV 2 (r)) , = W(r). (45) 

The above system of equations can be formally solved with the result: 



v 2 ( r) = -I/m dr <, 

r J r 

r 

V 1 (r) = -V 2 {r). 



(46) 
(47) 
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